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INTRODUCTION 


Precipitation,  clouds,  temperature,  humidity,  and  wind  affect  military 
tactics  significantly  in  many  areas,  e.g.,  performance  of  electro- 
optical  weapons  systems,  transport  and  diffusion  of  obscuration  mate- 
rials, toxic  and  incapacitating  chemical  agents  and  pollutants  in  the 
atmosphere,  and  maneuverability  of  aircraft. 

These  weather  variables  also  present  hydrologic  implications  to  military 
operations,  i.e.,  vehicular  trafficability,  to  the  extent,  or  degree, 
that  each  variable  influences  soil  moisture,  water  supply,  and  stream- 
flow.  At  least  90%  of  the  variation  in  water  runoff  from  a given  region 
is  caused  by  variations  in  the  intensity,  duration,  and  location  of 
rainstorms.  Snow  accumulation  is  a problem  because  it  is  so  difficult 
to  quantify;  but  even  if  snow  accumulation  is  known,  the  rate  of  melting 
assumes  major  importance  with  solar  insolation  becoming  a significant 
factor  probably  even  more  than  temperature.  Humidity  is  really  signi- 
ficant only  when  it  is  near  the  low  end  of  its  range  where  dehydration 
is  rapid  enough  to  be  troublesome.  Wind  velocity  coupled  with  low 
humidity  effectively  increases  the  rate  of  evapotranspiration. 

The  purpose  of  this  report  is  to  present  an  interpretative  review  of 
numerous  capabilities  for  measuring  and  forecasting  precipitation,  cloud 
cover  and  height,  and  profiles  of  temperature,  humidity,  and  windspeed 
that  exist  worldwide.  Emphasis  has  been  placed  upon  methodologies  and 
techniques  that  rely  on  the  collection  of  meteorological  data  by  remote 
means  and  upon  numerical  predictive  models  that  are  operational. 

There  has  been  no  intent  to  present  an  all-inclusive  treatment  of  the 
subject  matter.  Many  worthwhile  measurement  methods  and  forecasting 
techniques,  including  some  very  recent  developments,  have  not  been  cited 
due  to  the  limited  scope  of  this  study.  It  is  hoped,  however,  that  many 
omitted  methodologies  and  techniques  will  be  recognized  indirectly 
through  association  with  generalized  statements  or  similarity  to  more 
detailed  explanations  in  the  report. 

MEASUREMENT  CAPABILITIES 
Remote  Sensing  Methods 

Remote  sensing  methods,  i.e.,  those  involving  the  collection  of  data  by 
systems  not  in  direct  contact  with  the  phenomena  under  investigation, 
have  been  in  use  for  a number  of  years  in  the  form  of  weather  radar  (for 
detecting  precipitation  or  clouds)  and  other  techniques  [1].  Only 
recently,  though,  have  concerted  efforts  to  make  full  use  of  electro- 
magnetic spectrum  and  acoustic  probing  of  the  atmosphere  led  to  devel- 
opments in  Doppler  radar,  frequency-modulated,  continuous  wave  (FM-CW) 
radar  laser  systems,  and  acoustic  echo  sounders  that  are  helpful  in 
obtaining  data  on  all  the  important  parameters  of  the  lower  atmosphere 
[2]. 
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Remote  sensing  of  the  earth's  atmosphere  Is  often  hampered  by  the  pres- 
ence of  clouds  of  varying  density  over  the  region  being  surveyed  [3]. 
However,  by  the  end  of  this  decade.  It  Is  anticipated  that  satellite 
sounding  systems  will  be  able  to  penetrate  practically  all  cloud  condi- 
tions, exclusive  of  precipitating  clouds. 

Meteorological  Satellites 

Since  the  first  TIROS  satellite  was  launched  In  1960,  weather  predictions 
have  been  based  partially  on  Information  available  from  orbiting  weather 
satellites.  Originally,  artificial  satellites  were  used  simply  for 
pictorial  representations  of  cloud  distribution,  but  now  their  use  has 
extended,  via  the  remote  measurement  of  temperatures  of  land,  sea,  and 
cloud  surfaces,  to  the  reconstruction  of  complete  vertical  distributions 
of  temperature  and  humidity  [4]. 

In  addition  to  the  traditional  passive  sensors  operating  In  the  Infrared 
(IR)  and  visible  regions  of  the  electromagnetic  spectrum,  a phase  of 
satellite  sensors  has  now  begun  In  which  microwave  and  radar  systems  are 
being  practically  exploited  [5].  Radiometer  scanners  have  provided  ob- 
servations (with  Increased  Improvement  from  medium  to  very  high  resolu- 
tion) of  cloud  cover  and  severe  storms,  as  well  as  observations  from 
which  tropospheric  temperature  and  water  vapor  profiles  can  be  inferred 
C6]. 


US  Departrrent  of  Defense  Meteorological  Satellite  Program.  The  US 
Department  of  Defense  Meteorological  Satellite  Program  (DM§P)  Is  a 
unique  and  valuable  meteorological  data  system  [71.  The  sensors, 
communications,  and  data  processing  equipment  contribute  to  form  the 
most  responsive  operational  system  of  Its  kind.  Data  are  obtained  from 
space-borne  sensors  In  the  visual  and  IR  spectra.  Imagery  data  charac- 
teristics Include  both  visual  and  IR  sensor  resolution  capability  of 
0.33  nautical  mile  (nml)  for  limited  areas  and  2.0  nml  for  global  cover- 
age. The  spectral  bandwidth  of  the  visual  sensors  was  selected  to 
optimize  the  distinction  between  clouds,  ground,  and  water.  Electronic 
circuitry  In  the  sensors  converts  the  sensed  IR  energy  directly  to 
equivalent  blackbody  temperatures,  making  temperature  the  directly 
displayed  parameter.  The  sensitivity  of  the  2 nml  visual  channel  covers 
seven  orders  of  magnitude,  which  provides  useful  meteorological  infor- 
mation from  full  daylight  over  highly  reflective  scenes  to  an  Illumina- 
tion level  roughly  equivalent  to  half  moonlight.  One  additional  sensor 
furnishing  data  to  DMSP  Is  an  eight-channel  radiometer  for  the  deter- 
mination of  vertical  temperature  and  height  profiles.  Raw  data  are 
converted  Into  cloud  parameters  and  collated  with  conventional  meteoro- 
logical data  to  produce  a comprehensive  numerical  cloud  analysis.  With 
temperature  as  the  displayed  parameter,  the  resulting  enhancement  of  IR 
data  allows  the  pinpointing  of  storm  location  as  well  as  detailed  Infor- 
mation on  storm  cloud  structure. 

Currently,  DMSP  provides  data  sensed  near  dawn  and  near  local  noon  and 
their  nighttime  complements.  The  DMSP  communications  and  ground  pro- 
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cessing  system  are  designed  to  produce  a usable  product  within  5 minutes 
of  termination  of  the  data  stream.  For  direct  readouts,  a data  age  of  5 
to  20  minutes  when  ready  for  application  to  operational  decisions  is 
then  feasible.  Data  collected  by  the  satellite  are  both  broadcast  to 
military  receiving  terminals  in  the  satellite  local  area  and  stored  for 
later  broadcast  to  the  Air  Force  Global  Weather  Central  (AFGWC)  [8]. 

The  AFGWC  is  linked  to  its  readout  facilities  by  real-time,  wide-band 
communications. 

Soviet  Union  Meteorological  Space  System.  The  Soviet  Union  has 
operated  a meteorological  space  system,  ^'Meteor,"  since  the  late  1960's 
[9].  The  system  is  based  on  two  or,  more  often,  three  satellites  of 
similar  type  launched  into  almost  circular  polar  orbits  at  900  km 
altitude.  "Meteor"  satellites  carry  scanning  imaging  instruments  to 
obtain  cloud  images  and  snow  and  ice  cover  as  shown  by  temperature  con- 
trasts on  the  sunlit  and  dark  sides  of  the  earth,  as  well  as  to  measure 
the  radiation  temperature  of  clouds  and  earth  surfaces. 

The  Soviet  Union  plans  to  develop  a geostationary  satellite  system  over 
the  Indian  Ocean,  which  could  probably  be  launched  by  1978.  This  satel- 
lite is  expected  to  obtain  images  in  the  visible  and  IR  ranges,  collect 
information  from  instrumented  ground  platforms,  and  relay  data  between 
processing  centers. 


Weather  Radar 

As  a primary  information  instrument,  weather  radar  shows  the  shape  and 
dimensions  of  weather  disturbances  on  a screen.  Weather  radar  can  also 
provide  information  on  the  external  shape  and  physical  state  of  pluvial 
concentrations,  as  well  as  particle  volume  and  density  [10].  Several 
precipitation  types  (rain,  sleet,  snow)  have  been  distinguished  through 
variations  in  screen  brightness.  The  hydrometeor  detection  capabilities 
of  a radar  set  are  critically  dependent  upon  its  wavelength.  In  general, 
the  wavelength  required  to  detect  the  particles  decreases  with  decreased 
size  of  the  particles  [11]. 

Weather  Radar  Problems.  Weather  radar  application  can  be  quanti- 
tative,  as  well  as  qualitative,  in  the  measurement  of  precipitation; 
however,  quantitative  applications  require  correlation  between  precipi- 
tation and  instrument  response.  In  the  past,  a widespread  application 
of  weather  radar  to  hydrologic  studies  has  not  been  realized  because  the 
large  inherent  variance  of  the  radar  return  signal  from  rain  has  re- 
sulted in  uncertain  quantitative  measurements  of  rain  intensity  [12]. 

The  radar  presents  precipitation  intensity  distribution  in  a distorted 
form  because  of  attenuation  (reduction  in  intensity)  of  the  electro- 
magnetic (EM)  waves  in  precipitation  and  also  because  of  finite  dimen- 
sions of  the  sound  pulse  [13].  Cloud  and  rain  attenuation  has  to  be 
considered  at  wavelengths  below  the  10-cm  band  and  can  have  serious 
effects  in  the  3-cm  and  1-cm  bands. 
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Another  serious  obstacle  to  radar  estimation  of  rain  is  the  false  echo 
effect  resulting  from  anomalous  propagation  caused  by  density  variations 
in  the  atmosphere  [14]. 

In  1970,  important  initial  efforts  were  made  in  the  Soviet  Union  to 
improve  precipitation  and  cloud  measurements  by  means  of  either  two  sets 
of  radar  using  different  wavelengths  or  a single  set  of  radar  using  two 
wavelengths  to  determine  the  relative  absorption  sustained  by  the  two 
wavelengths  [151.  Numerous  investigations  have  been  made  by  different 
nations  during  the  past  few  years  relating  to  precipitation  intensity 
attenuation  and  depolarization  (due  to  nonsphericity  of  falling  rain- 
drops) at  microwave  and  millimeter  wavelengths. 

Much  effort  has  been  expended  to  find  paranieters  of  the  radar  echo 
distribution  that  are  usable  in  specifying  the  best  estimate  of  the 
radar  reflectivity-precipitation  rate  (Z-R)  relationship  for  each  storm 
[16].  There  have  been  numerous  studies  of  snow/EM  signal  interaction, 
but  they  predominantly  address  falling  snow  with  only  a few  addressing 
snow-covered  radar  targets  or  surfaces. 

Any  examination  of  the  capabilities  and  limitations  of  weather  radar  for 
quantitative  echo  intensity  measurements  must  consider  the  tradeoff 
among  beamwidth,  echo  integration  procedures,  and  scanning  time  limita- 
tions [17].  Fortunately,  new  alternatives  to  this  problem  have  been 
presented  by  recent  developments  in  electric  scanning  antennas,  fre- 
quency-agility noise-pulse  radars,  and  pulse  compression. 

US  National  Weather  Service  Weather  Radar  Program.  In  1976,  the  US 
National  Weather  Service  (NWS)  began  converting  their  weather  radars  to 
provide  an  automated  digital  data  stream.  Under  this  NWS  program, 
digital  video  integrators,  on-line  minicomputers,  and  necessary  auto- 
matic communications  equipment  will  be  added  to  existing  WSR-57  radars. 
Federal  Aviation  Administration  Air  Route  Traffic  Control  radars  are 
also  being  equipped  with  weather  and  fixed  map  units  which  provide  some 
weather  data  in  digital  forni. 

The  digital  video  technique  offers  advantages  for  the  digital  processing 
of  rainfall  data  for  flood  and  trafficabi 1 ity  forecasting.  It  is  also 
suited  for  producing  digitally  integrated  echo  intensity  contours  of 
severe  storms,  and  can  be  conveniently  transmitted  over  land  lines  for 
remote  display. 

Airborne  Weather  Radars.  Airborne  weather  radars  are  restricted  in 
size  of  antennas  carried  and  tend  to  use  shorter  wavelengths  that 
attenuate  appreciably  in  widespread  or  heavy  precipitation  [IS].  Unlike 
land-based  equipment  though,  these  radars  are  continually  moving  at  a 
high  rate,  resulting  in  opportunities  to  view  a given  piece  of  the 
precipitation  pattern  from  many  aspects. 

British  scientists  have  investigated  the  capability  of  airborne  side 
looking  radar  (SLR)  and  noted  the  following  system  deficiencies: 
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(1)  rain  and  snow  cause  system  attenuation;  (2)  SLR  resolution  is  not 
good;  (3)  imagery  is  small  scale,  and  distortion  and  shadowing  occur 
in  mountainous  areas  [19], 

Measurement  of  Precipitation 

Quantitative  measurements  of  precipitation  apparently  necessitate  the 
use  of  radar  having  a narrow  transmission  beam  in  all  planes  and  a 
wavelength  which  is  subject  to  negligible  attenuation.  A wavelength  of 
10-cm  seems  to  be  the  most  likely  to  reasonably  satisfy  all  of  these 
conditions,  although  relatively  large  antennas  and  relatively  high  peak 
transmitted  power  are  then  required.  At  wavelengths  less  than  10-cm,  no 
precise  measurement  can  be  made  without  a correction  being  made  to  the 
power  returned  to  the  radar.  As  a compromise,  C-band  (5.5-cm  wave- 
length) radars  offer  the  advantage  of  nominal  precipitation  attenuation 
combined  with  good  detection  capability  for  frozen  and  low  rates  of 
precipitation  plus  adequate  beamwidth  resolution  of  the  echoes. 

Methods  of  Estimating  Rainfall 

a.  Satellite  Data 

Many  methods  have  concentrated  on  estimating  convective  rainfall.  In 
these  estimation  methods,  regions  of  convective  activity  must  be  iso- 
lated from  background  cloudiness  and  the  extent  and  significance  of 
nonconvective  rainfall  must  be  determined.  Methods  using  visible 
satellite  data  have  tackled  the  first  difficulty  through  an  enhancement 
technique  which  suppresses  regions  of  shallow  convection  and  inactive 
stratiform  clouds  by  setting  a cloud  brightness  threshold.  Typically, 
radiance  response  on  meteorological  satellites  is  set  such  that  the 
upper  10%  to  40%  of  the  brightness  range  contains  the  regions  of  deep 
convection.  Similarly,  IR  techniques  try  to  discern  regions  of  con- 
vective activity  and  estimates  of  rainfall  by  comparing  equivalent 
blackbody  temperatures  to  climatological,  regionally  arranged,  or  actual 
temperature  soundings.  Both  methods  are  flawed.  Brightness  enhancement 
suffers  from  signal  saturation  and  dependence  on  sun  angle  and  viewing 
geometry;  while  the  IR  method  is  compromised  by  cirrus  contamination  and 
resolution. 

By  comparing  satellite  brightness  contoured  cloud  masses  with  contoured 
radar  echoes,  it  appears  possible  to  use  satellite  pictures  for  infer- 
ring precipitation  frequency  (occurrence),  extent,  intensity,  and  rate 
of  change. 

A method  developed  in  the  Soviet  Union  visibly  identifies  precipitation 
and  gives  an  approximate  estimation  of  its  intensity  using  satellite 
microwave  measurements  on  two  orthogonal  polarizations  at  the  0.8-cm 
region  [20].  Developed  from  computational  data  of  microwave  radiation 
intensity,  this  method  is  based  on  the  dependence  of  absolute  values  of 
radio  brightness  temperature  measured  on  vertical  and  horizontal  polar- 
izations and  their  difference  on  total  content  of  water  vapor  and  liquid 
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water  in  the  atmosphere.  This  method  makes  it  possible  to  identify 
cloud  formations  and  zones  of  precipitation,  and  to  approximate  the 
integral  water  content  of  clouds  and  rainfall  intensity. 

b.  Other  Remote  Techniques 

Dual-  and  multi -Doppler  radars  have  been  used  to  specify  the  dimensional 
particle  velocities  within  convective  storms  as  related  to  heavy  rain- 
fall. The  drop  size  distribution  in  rain  (with  difficulties  involved  in 
the  case  of  snow)  can  be  computed  from  the  Doppler  spectrum,  provided 
the  updraft  velocity  can  be  estimated  [211.  Definite  seismic  and 
acoustic  intrusion  sensor  activation  patterns  have  also  been  correlated 
with  various  rainfall  rates  [223. 

Monostatic  lidar  has  been  explored  as  a means  of  determining  the  rain- 
fall rate  over  an  extended  atmospheric  path  [23].  In  accounting  for  the 
multiple  scattering  of  light  in  rain,  estimates  have  been  derived  for 
the  optical  extinction  coefficient  of  rain  as  a function  of  range  from 
the  lidar  returns.  Lidar  soundings  of  the  atmosphere  under  a variety  of 
conditions  (snowfall,  rain,  drizzle,  etc.)  have  given  sharp  polarization 
differences  of  signals  that  can  be  used  to  identify  the  meteorological 
characteristics  of  hydrometeors  [24]. 

The  Environmental  Research  Laboratory,  N0W\,  has  been  probing  various 
weather  systems  by  use  of  multiwavelength  lidar  combined  with  shortwave 
radar  and  IR  radiometry  [25].  Of  particular  interest  is  a laser  probing 
system  that  uses  a flashlamp  pumped  dye  laser  as  the  radiation  source 
and  a plastic  (Fresnel)  lens  in  the  receiver  to  collect  the  radiation 
backscattered  by  the  atmosphere  [26].  This  system  is  claimed  to  be  less 
expensive,  more  reliable,  and  more  compact  than  systems  incorporating 
ruby  lasers  and  more  conventional  receiver  optics. 

c.  Rain  Gages 

The  tipping  bucket  precipitation  gage,  weighing- type  rain  gage,  and  the 
Jardi  rate-of-rainfal 1 gage  are  just  a few  of  the  many  precipitation 
measuring  gages  used  throughout  the  world.  However,  during  recent 
years,  new  improved  instrumentation  has  been  developed  to  increase  the 
accuracy  of  precipitation  measurements.  The  Workman  intensity  gage 
which  senses  the  flow  of  water  down  a movable  trough  and  the  Raymond- 
Wilson  gage  which  measures  the  electrical  resistance  of  flowing  water 
are  both  superior  to  standard  gages  because  of  their  rapid  response  time 
and  capacity  to  measure  high  intensity  showers  continuously  and  auto- 
matically [27J.  A recording  rain  gage,  with  electronic  sensors  and 
circuitry  to  provide  a sensitivity  of  0.1  mm  in  recording,  permits  the 
accurate  sensing  of  mists  and  drizzle,  as  well  as  rainfall  of  almost  any 
known  rate  [28].  A laser  rain  gage  measures  precipitation  by  the 
scattering  of  light  from  raindrops.  The  sampling  medium  is  a collimated 
beam  from  a helium-neon  laser.  The  amount  of  light  scattered  is  a 
function  of  the  number  and  size  of  drops  intercepting  the  beam;  an 
equation  relating  rainfall  rate  and  scattered  light  is  derived. 
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Some  scientists  have  visualized  that  radar  (which  can  observe  precipita- 
tion over  a very  large  area)  ultimately  will  replace  rain  gages.  Others 
have  taken  the  position  that  radars  should  complement  rain  gages  and 
that,  when  used  together,  they  will  provide  significantly  better  infor- 
mation on  rainfall  than  can  be  obtained  from  a rain  gage  network  alone. 

In  general,  investigators  have  shown  that  radars,  properly  calibrated 
with  a rain  gage  in  the  area,  can  make  reasonably  accurate  measurements 
of  rainfall  intensity  and  accumulated  rainfall  over  a small  watershed. 

Snow  Measurements.  A disdrometer  technique  for  measuring  snowflake 
size  distributions  uses  a pulsed  ruby  laser  which  stops  particle  motion 
and  provides  instantaneous  periodic  samplings  over  the  history  of  the 
particles  [29J.  An  acoustic  device  that  tallies  all  particles  over  20ym 
provides  an  accurate  count  of  large  snowflakes  falling  into  the  sensor 
1301. 

A new  type  of  dilatometer  has  been  designed  to  determine  the  free  water 
content  of  wet  snow  within  an  error  of  1%  from  the  decrease  in  volume 
caused  by  the  process  of  melting  snow  [311. 

The  use  of  natural  gamma  radiation  from  the  soil  as  a basis  for  snow- 
water equivalent  measurements  at  remote  sites  has  been  under  investiga- 
tion by  the  NWS  since  1970  [321.  Results  have  indicated  measurements 
with  about  5°'  error  in  the  5-40  cn,  water  equivalent  range  during  periods 
uncomplicated  by  precipitation  or  considerable  change  in  soil  moisture. 

It  has  also  been  found  that  airborne  measurement  of  the  areal  average  of 
the  water  equivalent  of  snow  by  the  terrestrial  gamma  attenuation  method 
is  accurate  to  1.2-cm,  or  10%,  whichever  is  greater  [33].  This  method 
has  good  potential  for  large,  rapid,  areal  measurements  of  water  equiva- 
lent over  flat  or  rolling  terrain.  The  attenuation  of  highly  penetrating 
cosmic  radiation  appears  to  have  excellent  potential  for  measuring  the 
water  equivalent  of  snow  cover  in  extremely  deep  snow.  A two-scintilla- 
tion detector  setup,  one  above  and  one  beneath  the  snow,  has  produced  a 
water  equivalent  measurement  accuracy  of  better-  than  1%  in  measuring 
100  cm  of  water  with  a 24-hour  measurement  time  [34]. 

Remote  sensing  of  snow  cover  by  SLR,  with  wavelengths  of  1 cm  or  smaller, 
has  provided  quick,  comprehensive  information  about  snow  depth,  density, 
and  water  content  [35]. 

Visible  and  near-visible  light  sensors  carried  on  meteorological  satel- 
lites can  also  measure  the  snow  area  extent,  but  cannot  routinely  detect 
snow  under  clouds  or  a forest  canopy.  Furthermore,  much  research  needs 
to  be  done  on  the  electrical  (including  scattering)  properties  of  snow 
before  efficient  all-weather  remote  sensing  systems  can  be  designed  [36]. 

Hail  Sensing  Instruments.  Surface  in  situ  hail  sensing  instruments 
have  been  classified  into  two  major  categories  [37].  The  integrating 
sensors  include  the  hailpad,  hailstool,  hail  cube,  and  hail-wind  detector. 
Time-recording  sensors  include  the  Geophone  hailgage,  Illinois  Water 
Survey  hailgage,  HOAA  hail  momentum  sensor,  and  a hail-rain  separator. 
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The  integrating  sensors  are  relatively  inexpensive  and  yield  useful 
data,  including  estimates  of  wind  effects.  The  recording  type  is  10-100 
times  more  expensive,  but  for  this  additional  cost,  obtains  detailed 
time  information  that  is  not  available  from  integrating  instruments. 

A recent  developmental  effort  in  remote  hail  detection  has  been  the 
dual -wavelength  radar  hail  detection  technique  proposed  by  Eccles  and 
Atlas  [37].  Based  on  the  variable  scattering  properties  of  hail  at 
different  wavelengths,  this  approach  has  met  with  limited  success. 

Another  radar  hail  detection  technique  is  based  on  the  depolarization  of 
microwave  radiation  by  nonspherical  hailstones.  Tests  of  this  approach 
in  Canada  have  lent  support  to  its  validity. 

A theoretical  investigation  has  shown  that  the  ratio  of  radar  reflectiv- 
ity factors  at  10  cm  and  3 cm  wavelengths  is  useful  in  delineating  hail 
regions  and  can  also  be  used  to  estimate  mean  hail  diameter  under 
certain  conditions  [38].  Regions  of  hail  can  be  identified  as  regions 
of  large  ratios  after  taking  attenuation  into  account  where  necessary. 

Measurement  of  Cloud  Amount  and  Height 

Weather  Satellite  Applications.  The  application  of  weather  satel- 
lites  has  been  a milestone  in  the  monitoring  of  large-scale  clcuciness, 
even  though  the  utility  of  meteorological  satellite  data  is  currently 
limited  because  of  an  inability  to  precisely  define  the  heights  and 
small-scale  structure  of  clouds.  Satellite  cloud  photographs  have  been 
used  to  determine  the  approximate  altitude  of  clouds  by  applying  an 
approximate  equation  that  relates  cloud  altitude,  angle  of  sun  elevation, 
and  cloud  image  to  corresponding  shadow  distance  as  measured  directly  on 
the  photograph  [39].  Another  rapid  objective  method  calculates  the 
horizontal  dimensions  of  cloud  systems  based  upon  the  height  of  satel- 
lite orbit,  dimensions  of  the  picture  taken,  and  view  angle  of  the  TV 
camera  [40]. 

Satellite  digital  images  have  also  been  used  to  study  the  time-changing 
cloud  and  thermodynamic  structure  of  cloud  clusters  and  their  environ- 
ment over  the  oceans  [41], 

Characteristic  features  and  temporal  variations  of  clouds  have  been 
obtained  from  cloud  brightness  data  and  radar  plan-position  indicator 
(PPI)  pictures.  From  the  PPI  display,  one  can  immediately  read  the 
range  and  bearing  to  the  target.  A comparison  of  radar  heights  of  large 
cumulus  clouds  to  their  reflected  solar  radiance  has  shown  strong  cor-" 
relation  (0.88),  thus  suggesting  that  cumulus  cloud  heights  and  growth 
may  be  inferred  from  a satellite  platform  without  use  of  ground-based 
radar  [42]. 

Other  Methods  for  Evaluating  Cloud  Amount.  An  automatic  cloud- 
amount  recorder  developed  by  the  Swedish  Kieteorological  and  Hydrological 
Institute  is  based  on  the  principle  that  each  cloud  elem.ent  is  charac- 
terized by  a certain  temperature  at  the  cloud  base  and  a certain  aper- 
ture angle  from  the  place  of  observation  143].  By  measuring  the  tempera- 
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ture  and  knowing  the  temperature  profile  from  the  latest  radiosonde 
measurements,  it  is  possible  to  decide  whether  the  cloud  is  a high, 
middle,  or  low  cloud  and  to  determine  the  cloud  size.  The  model  re- 
corder, equipped  with  rotating  aluminum  mirrors  and  IR  thermometer, 
scans  from  90  to  40  degress  above  the  horizon.  Three  different  temper- 
ature layers  (corresponding  to  high,  middle,  and  low  cloud)  are  selected 
at  three  potentiometers,  giving  three  reference  voltages.  The  signal 
from  the  IR  sensor  is  fed  to  a level  detector  with  three  outputs,  which 
indicate  that  the  cloud  is  detected  at  the  corresponding  cloud  level. 

Out  of  such  signals,  the  amounts  (in  eights  of  sky  coverage)  of  high, 
middle,  and  low  cloud  are  integrated.  Reasonably  good  results  were 
obtained  as  compared  with  manual  observations. 

Other  instruments  have  also  been  built  to  evaluate  cloud  amount  C44]. 

One  is  based  on  the  difference  in  the  IR  emission  of  clouds  and  clear 
sky  and  uses,  as  a sensor,  a net  radiometer  with  screening  consisting  of 
a thin  metal  sheet  blackened  on  the  internal  face.  Another  exploits  the 
intensity  of  sky  light  in  the  two  extreme  visible  regions.  Two  filters 
and  two  sensors  are  used  to  evaluate  the  fraction  of  polarized  light  in 
sky  light.  A sky  mirror  system  has  been  used  to  take  pictures  of  the 
whole  sky,  with  the  clouds  being  projected  perpendicularly  upon  the 
plane  of  the  horizon  to  form  a cloud  cover  map  [453.  All -sky  camera 
photographs  have  been  used  to  determine  low,  middle,  high,  and  total 
cloud  cover,  in  addition  to  diurnal  cloud  variations  over  ocean  areas. 

Measurement  of  Cloud  Height.  Various  devices  have  been  used  for 
determining  cloud  height  L46JT  fiany  stations  use  automatic  ground-based 
ceilometers,  e.g.,  fixed  beam  ceilometer,  rotating  beam  ceilometer  and 
an  eyesafe  laser  ceilometer,  to  make  cloud  height  measurements;  with 
manual  techniques,  e.g.,  human  estimation,  ceiling  light,  and  ceiling 
balloons  being  used  as  secondary  systems  [473.  Cloud  height  measuring 
equipment  at  the  Berlin  (Germany)  Meteorological  Institute  includes  an 
X-band  vertical  radar  (former  PPI  radar),  a stereotelemeter,  and  a 
ceilograph  (ceiling  projector  with  pulsed  light)  [48]. 

In  1972,  the  Japanese  announced  the  testing  of  a prototype  laser  radar 
system  which  transmits  a pulsed  laser  light  and  emits  a much  shorter 
wavelength  than  that  of  the  conventional  microwave  system  [49].  It  was 
claimed  that  this  system  is  capable  of  producing  scope  images  of  both 
thick  and  thin  clouds,  cloud  thickness,  fog,  rain  and  snow.  As  related 
to  laser  investigation  of  clouds,  the  influence  of  molecular  and  ^rosol 
scattering  and  of  molecular  absorption  upon  the  magnitude  of  reflected 
echo  signals  has  been  evaluated  [503. 

The  US  Air  Force  Geophysics  Laboratory  has  evaluated  a ruby  lidar  ceil- 
ometer for  applicability  to  cloud  heights  [513.  It  was  concluded  that 
lidar  indicates  an  accurate  presentation  of  cloud  structure  and  is 
superior  to  a standard  rotating  beam  ceilometer  as  a cloud  height 
measuring  device.  The  US  Army  has  been  developing  a portable  cloud 
ceiling  and  visibility  sensor  through  application  of  the  AN/GVS-5  Laser 
Range  Mnder.  Cloud  heights  are  measured  by  reference  to  the  time  lapse 
from  transmitted  to  detected  signal  [523. 
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Airborne  IR  scanners  have  also  been  used  to  obtain  stereo  imagery  of 
clouds  from  which  cloud  forms  and  stratification,  including  cloud  base 
and  top  heights,  have  been  determined. 

Measurement  of  Temperature,  Humidity,  and  Windspeed  Profiles 

Temperature  Profiles 

a.  Radiance  Measurements 

Ik  radiance  measurements  to  obtain  absorption  and  emission  data  for 
temperature  and  humidity  profiles  have  been  made  by  sensors  based  on  the 
ground  and  water  (buoys)  as  well  as  sensors  carried  aloft  by  aircraft, 
balloons,  high-altitude  rockets  and  satellites  [531.  Vertical  tempera- 
ture profiles  are  derived  from  satellite  radiation  measurements  by 
inverting  the  integral  form  of  the  radiative  transfer  equation  [541. 
Instrumentation  for  IR  soundings  on  the  basis  of  thermal  radiation 
includes  the  Michelson  interferometer,  gravity  spectrometers,  filter 
radiometers,  and  microwave  detectors  [551.  In  the  limb  scanning  tech- 
nique, emitted  radiation  from  the  atmosphere  is  measured  by  a radiometer 
as  it  scans  across  the  planetary  limb,  with  the  measurements  being 
inverted  to  obtain  the  atmospheric  temperature  structure  [561. 

The  solar  occultation  technique,  based  upon  measurements  of  the  sun  tan- 
gentially through  the  earth's  atmosphere  from  a satellite,  provides 
profiles  of  upper  air  temperature  (by  inversion)  from  absorption  spectra 
data  for  niinor  atmospheric  constituents  [571.  Vertical  resolution  and 
sensitivity  have  been  found  to  be  greater  for  the  occultation  technique 
than  for  the  radiance  technique  of  measurement. 

A ground-based  IPv  spectroradiometer  has  been  used  to  measure  the  verti- 
cal temperature  profile  of  the  lower  atmosphere  from  the  surface 
to  6 km.  The  resultant  profile  had  an  accuracy  comparable  to  that  of 
radiosondes  !58]. 

Temperature  profiles  have  also  been  inferred  from  ground-based  radio- 
metric  observations  of  the  emission  by  atmospheric  oxygen  through  a 
quasi -linear  integral  equation  relationship  [591. 

b.  Lasor-Lidar  Methods 

Measurements  of  atmospheric  temperature  profiles  have  used,  as  a physical 
basis  for  the  measurements,  the  temperature  sensitivity  of  the  envelope 
of  the  Raman  rotational  hand  backscattered  from  the  N2  and  O2  content 
[60].  The  measurements  were  made  by  a laser  radar  which  used  a dual 
optical  channel  configuration  in  the  receiver;  thus,  the  backscattered 
radiation  is  split  between  the  two  channels,  each  of  which  contains  a 
narrow  band  optical  filter  designed  to  pass  a prearranged  portion  of  the 
rotational  spectrum.  A typical  lidar  differenced  signal  (scope  trace) 
is  presented.  Future  plans  call  for  implementing  the  lidar  with  a 
transient  recorder  - minicomputer  receiver  package  to  provide  not  only 
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the  sensitivity  needed,  but  also  on-line  temperature  profiles. 

A large  laser  radar  system  has  been  built  for  use  in  measuring  Rayleigh 
and  Raman  scattering  in  two  colors  from  the  atmosphere  [611.  The  two 
color  capabilities  of  the  system  are  used  to  obtain  data  on  the  differ- 
ences between  Rayleigh  and  Mie  scattering  which  exhibit  different 
scattering  functions  of  different  wavelengths.  This  laser  radar  has 
proven  to  be  effective  in  obtaining  data  on  clouds,  layering  of  water 
vapor,  and  other  inhomogeneities  that  are  difficult  to  measure. 

The  thermodynamic  distribution  within  a molecular  absorption  band  of  a 
gas  has  been  probed  by  means  of  lidar,  although  relatively  weak  absorp- 
tions are  involved  [621.  Another  laser  method  of  determining  the  atmo- 
spheric temperature  distribution  is  based  on  measurements  of  the  combina- 
tional light  scattering  ratio  of  Stokes  to  anti-Stokes  constituents  [631. 

c.  Acoustic  Techniques 

Spectral  broadening  of  acoustic  and  radio  waves  scattered  by  atmospheric 
turbulence  has  been  studied,  with  special  emphasis  on  Doppler  radar 
applications  [641.  A single  formalism  has  been  derived  for  the  scattering 
of  both  acoustic  and  electromagnetic  vvaves  which  accounts  for  the  spectral 
broadening  resulting  from  the  deformation  of  the  turbulent  medium,  the 
finite  size  of  the  observed  volume,  the  bulk  motion  of  air  through  this 
volume,  and  the  pulse  modulation.  The  backscatter  spectrum  associated 
with  each  of  these  broadening  mechanisms  has  also  been  computed. 

A relatively  new  technique  for  the  remote  measurement  of  temperature  in 
the  lower  atmosphere  employs  the  Radio  Acoustic  Sounding  System  (RASS) 
[651.  A burst  of  sound  propagating  upward  in  the  air  is  tracked  by 
Doppler  radar.  Air  temperature  at  each  height  is  determined  from  the 
instantaneous  speed  of  the  sound  pulse,  and  a complete  profile  can  be 
obtained  in  a few  seconds.  The  RASS  penetrates  clouds,  heavy  pollution, 
precipitation,  and  low-level  inversions.  It  has  a temperature  error 
much  less  than  1®C  from  near  ground  level  to  1-3  km  for  most  wind  condi- 
tions, and  is  a much  more  cost-effective  procedure  than  the  radiosonde 
technique. 

Backscatter  echoes  from  acoustic  tone  bursts  have  been  compared  with 
temperature  profiles  measured  by  radiosondes  through  the  height  range  of 
50-700  m [661.  In  tests,  the  sounder  recorded  a somewhat  lower  tempera- 
ture inversion  than  radiosonde  data  indicated,  probably  because  the 
sounder  has  an  advantage  of  continuously  monitoring  the  inversion  struc- 
ture and  undulations  produced  by  wind  shear. 

Humidity  Profiles.  The  measurement  of  atmospheric  water  vapor 
using  the  Raman  backscatter  from  a pulsed  laser  is  feasible  [67];  and 
the  possibility  of  deriving  humidity  profiles  from  the  backscatter  of 
acoustic  radar  signals  has  been  confirmed  [68].  Humidity  profiles,  have 
also  been  determined  by  lidar  measurements  within,  and  adjacent  to,  a 
water  vapor  absorption  line  [691. 
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But  consideration  of  the  effects  of  temperature  and  pressure  on  the 
determination  of  water  vapor  densities  from  ruby  laser  radar  returns  has 
brought  forth  the  suggestion  that  water  vapor  densities  be  deduced  from 
the  line  strength,  which  is  independent  of  pressure,  rather  than  from 
the  absorption  coefficient  [70].  However,  the  accuracy  of  measuring 
water  vapor  concentration  and  atmospheric  humidity  by  laser  radiation 
scattering  is  dependent  upon  the  errors  of  measuring  the  scattered 
radiation  intensity  and  radiation  wavelength  [711.  The  inference  of 
water  vapor  profiles  from  ground-based  remote  measurements  of  spectral 
radiances  has  indicated  that  a fine-scale  detailed  profile  structure 
cannot  be  reconstructed  through  the  linear  inversion  technique  which 
linearizes  the  radiative  transfer  function,  but  a nonlinear  iteration 
technique,  which  minimizes  the  root  mean  square  residual  of  the  random 
noise  along  the  direction  of  steepest  descent,  can  be  used  to  infer  a 
reasonably  stable  solution  [72  1. 

Wind  Profiles 


a.  Loran  NAVAID  System 

Rawinsondes,  dropsondes,  and  pilot  balloon  (pibal)  observations  have  all 
been  used  routinely  for  many  years  to  obtain  vertical  wind  profiles. 

But  to  provide  more  accurate  wind  measurements,  equipment  of  more  recent 
vintage  is  available.  The  Loran  flAVAID  system  for  measuring  the  vertical 
profile  of  vnnd  is  capable  of  a 1.4-knot  accuracy,  and  the  Onega  equip- 
ment an  accuracy  of  4.3  knots,  within  a 600-m  height  interval  (731. 

b.  Dropwindsonde  System 

During  the  GARP  Atlantic  Tropical  Experiment  (GATE)  in  1974,  a unique 
windfinding  system,  called  a dropwindsonde,  was  used  in  field  operations 
for  the  first  time  [741.  Deployed  from  an  aircraft  flying  at  about 
10  km  altitude,  the  dropwindsonde  transmitted  vertical  profiles  of 
horizontal  winds,  pressure,  temperature,  and  humidity,  while  descending 
on  a parachute  at  the  rate  of  25  mb/niin.  The  sonde  has  four  sensors, 
including  a bead  thermistor  for  temperature,  a carbonhygristor  for 
relative  humidity,  and  an  improved  aneroid  cell  for  pressure.  Wind- 
finding is  accomplished  by  phase  measurements  of  Omega  navigation 
signals  received  and  retransmitted  by  the  sonde  to  the  aircraft. 

c.  US  National  Weather  Service  "NEXAIR"  System 

During  the  past  few  years,  the  NWS  has  made  considerable  progress  in 
developing  a new  upper  air  sounding  system  (NEXAIR)  that  will  satisfy 
projected  requirements  for  large-  and  small-scale  upper  air  data  [75], 
Using  navigational  aids  and  a nonrotating  antenna  system  under  mini- 
computer control,  403-MHz  frangible  sondes  will  relay  pressure,  tempera- 
ture, and  humidity  data  to  the  ground  by  telemeter.  NEXAIR  promises  to 
be  a safe,  cost-effective,  low-maintenance  system  capable  of  being 
operated  in  a wide  variety  of  extreme  environments. 


d.  Jimsphere  Technique 


The  Jimsphere  technique,  which  involves  the  tracking  of  a spherical 
balloon  by  a high  precision  radar  (FPS-16)  system,  has  provided  an 
improved  method  for  obtaining  vertical  wind  profile  data  at  25-m  inter- 
vals from  near  the  ground  to  18  km  [76].  As  an  imorovement  to  this 
system,  a lightweight  temperature  sensor  (Jimsonde)  was  developed  to 
give  high  resolution  temperature  data  for  each  25-m  height  interval 
simultaneously  with  the  wind  data  [77]. 

To  measure  wind  profiles  more  rapidly,  and  to  improve  spatial  and  tem- 
poral resolution,  new  devices  are  being  developed  that  use  electromag- 
netic and  acoustic  waves  to  measure  winds  remotely. 

e.  Doppler  Techniques 

Vertical  wind  can  be  measured  acoustically  with  a single  axis  monostatic 
system.  The  Doppler  shift  of  the  returned  signal  is  related  to  the 
vertical  velocity  of  the  scattering  volume,  which  is  the  vertical  com- 
ponent of  the  wind.  By  tilting  the  single  monostatic  antenna,  a com- 
ponent of  the  horizontal  wind  can  also  be  measured.  A weakness  in 
making  measurements  with  a monostatic  system  is  that  the  return  signals 
are  not  always  continuous  under  some  types  of  thermal  structure.  The 
use  of  a bi static  system,  where  transmitter  and  receiver  are  separated, 
helps  to  solve  this  problem. 

The  acoustic-Doppler  technique  measures  the  motion  of  wind-carried 
eddies  from  the  Doppler  shift  of  backscattered  acoustic  pulses.  It  can 
measure  wind  profiles  up  to  0.5  km,  but  is  ineffective  in  precipitation 
and  in  high  winds.  Acoustic  absorption  also  imposes  height  limits. 

Doppler  lidars  measure  the  velocity  of  wind-advected  aerosols  and  so 
depend  on  the  presence  of  such  aerosols.  Furthermore,  they  are  vulner- 
able to  attenuation  by  precipitation.  Another  problem  is  that  CW  Doppler 
lidars  have  a range  resolution  that  degrades  as  the  range  is  increased. 

The  METRAC  Positioning  System  is  a self-containec,  ground-based,  radio- 
location system  that  uses  the  Doppler  principle  to  track  a lightweight, 
expendable  transmitter.  The  purpose  of  the  system  is  to  provide  an 
economical,  versatile,  and  accurate  technique  for  tracking  meteorological 
balloons  [78].  In  field  tests,  by  tracking  the  same  balloons  simultan- 
eously with  a prototype  METRAC  system,  the  rawinsonde  system,  and 
theodolite  equipment,  it  was  concluded  that  the  METRAC  system  is  capable 
of  providing  high  resolution  wind  profiles. 

A Doppler  radiosonde  system  called  the  Safesonde  is  a minimum  weight, 
mininum  cost,  and  minimum  hazard  sounding  system  for  measuring  the  winds 
from  the  surface  to  15  km  [79].  This  system  consists  of  a simple  balloon- 
borne  transmitter  tracked  by  five  receivers  to  provide  position  data. 


f.  FM-CW  Microwave  Radar  Technique 

Another  new  technique  for  measuring  winds  in  the  lower  atmosphere  with  a 
microwave  radar  is  an  extension  of  the  well  known  FH-CW  technique,  and 
it  has  the  same  advantages  of  relatively  low  cost  and  high  flexibility 
for  a clear-air  radar.  The  microwave  technique,  in  principle,  is  immune 
to  most  of  the  factors  that  limit  the  usefulness  of  the  optical  and 
acoustic  techniques.  This  technique  offers  an  all-weather  measurement 
capability;  ;inds  in  the  clear  air  can  be  measured  to  heights  in  excess 
of  1 km,  using  microwave  FM-CW  radar  (10-cm  wavelength  and  two  identical 
8-foot  parabolic  antennas  for  transmitting  and  receiving  the  signal) 
with  Doppler  processing  [80].  With  higher  power  and  larger  antennas, 
there  should  be  a capability  to  probe  greater  heights.  Aside  from 
system  parameters,  the  maximum  radar  range  depends  on  the  availability 
of  turbulent  scatterers.  Nevertheless,  radar  wind  measurement  offers 
significant  improvements  in  spatial  and  temporal  coverage  over  in  situ 
techniques. 

Ground-Based  Remote  Sensing  Limitations 

The  applicability  of  multi-frequency  backscatter  acoustic  radar  for  moni- 
toring vertical  humidity  profiles,  Doppler  radar  for  measuring  temperature 
and  wind  velocity,  and  polarized  bistatic  backscatter  lidar  for  measuring 
temperature  and  air  density  is  currently  under  investigation. 

However,  it  should  be  pointed  out  that  important  limitations  currently 
exist  in  the  ground-based  remote  sensing  of  meteorologically  significant 
parameters  [81].  For  example,  most  optical  and  IR  measurement  systems 
have  greatly  reduced  measurement  capabilities  in  fog,  cloud,  or  precipi- 
tation. On  the  other  hand,  some  radar  systems  require  the  presence  of 
precipitation  (or  chaff)  to  provide  echoing  targets.  Acoustic  systems 
are  usually  adversely  affected  by  the  impact  noise  produced  by  rain  or 
hail,  but  (unlike  optical  systems)  they  are  not  adversely  affected  by 
cloud  or  fog. 


FORECASTING  CAPABILITIES 

Two  differing  viewpoints  exist  concerning  the  laws  of  atmospheric  move- 
ments and  how  to  forecast  them.  One  viewpoint  is  based  on  the  assump- 
tion that  the  laws  of  atmospheric  movement  are  relatively  fixed  and  that 
with  enough  observational  data,  a set  of  equations  on  atmospheric  motion 
could  be  computerized  and  accurate  weather  forecasts  produced  - this  is 
numerical  forecasting.  The  other  viewpoint  entails  the  premise  that  the 
laws  of  atmospheric  motion  are  random,  and  that  regardless  of  the  size 
of  the  observational  base,  atmospheric  processes  cannot  be  completely 
described  and  probability  methods  must  be  used  - this  is  probability 
forecasting. 


Numerical  Prediction  Models 

The  simplest  atmospheric  models  used  in  numerical  weather  prediction  are 
based  on  the  assumption  of  a barotropic  atmosphere  [82].  Individual 
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barotropic  models  differ  with  regard  to  structure,  but  have  substantially 
the  same  characteristics  in  that  they  do  not  include  the  transformation 
of  potential  into  kinetic  energy. 

Most  numerical  models,  for  extended  range  prediction,  use  primitive 
equations  of  hydrothermodynamics  that  are  integrated  over  either  the 
hemispheric  or  global  domain  with  varying  vertical  resolution  [83]. 
Depending  on  their  purpose  and  current  state  of  development,  the  models 
differ  considerably  in  their  way  of  treating  the  physical  processes. 
Different  finite-difference  approximation  techniques  are  used  with 
various  degrees  of  parameterization  of  the  atmospheric  boundary  layer, 
cumulus  convection,  radiation  and  topography.  Recent  advances  have 
included  the  development  of  numerical  models  suitable  for  studying  the 
general  atmospheric  circulation  for  application  to  extended-range  pre- 
dictions. 

Semi-implicit  schemes  (for  increasing  the  effectiveness  of  computational 
forecasts  for  5-7  days  and  more)  have  also  been  applied  to  the  integra- 
tion of  a system  of  primitive  equations  describing  a baroclinic  model 
which  recognizes  atmospheric  energy  conversion  from  potential  to  kinetic 
energy  [84]. 

US  National  Meteorological  Center  Models.  Since  1966,  a hemispheric 
six-layer  baroclinic  model  which  uses  the  primitive  (Newtonian)  equations 
(PE)  has  been  used  as  the  National  Meteorological  Center  (NMC)  operational 
model  [85],  A regional  Limited-Area  Fine-Mesh  (LFM)  model  was  added  in 
1971,  The  forecasting  community  in  the  US  now  receives  forecasts  out  to 
48  hours  from  the  barotropic  hemispheric  PE  and  regional  LFM  models  twice 
each  day,  based  on  analyses  made  from  0000  and  1200  GMT  observations.  The 
PE  model  is  also  run  out  to  84  hours,  once  per  day,  from  the  0000  GMT 
observations. 


Use  of  Statistics  and  Cl iinatology  in  Forecasting 

Numerical  prediction  is  just  one  weapon  in  the  forecaster's  arsenal; 
statistical  methods,  extrapolation  techniques,  and  empirical  rules  will 
continue  to  have  their  place  in  weather  forecasting  for  the  foreseeable 
future. 

In  recent  years,  synoptic-statistical  and  hydrodynamical  methods,  along 
with  statistical  climatology,  have  been  applied  to  both  short-range  and 
long-range  forecasting.  Modern  developments  in  statistical  meteorology 
have  included  the  areas  of  stochastic-dynamic  prediction,  statistical 
weather  forecasting,  probability  forecasting  and  time  series  analysis 
[86].  Short-range  forecasting  approaches  have  been  based  on  statistics, 
local  climatology,  a combination  of  climatology  and  persistence,  and  the 
use  of  multiple  prediction  statistics,  i.e.,  stepwise  regression  and 
forecasting  from  statistically  determined  graphs. 
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The  rationale  of  a stochastic-dynamic  approach  to  numerical  weather 
prediction  deals  with  an  initial  uncertainty  by  considering  an  infinite 
ensemble  of  initial  states  in  space-phase,  with  relative  frequencies 
within  the  ensemble  being  proportional  to  probability  densities  [87]. 

The  evolution  of  the  ensemble  in  time,  as  given  by  the  stochastic- 
dynamic  equation  set,  is  based  upon  the  original  deterministic  hydro- 
dynamic  equation  set. 

Climatology  has  always  formed  an  integral  part  of  weather  forecasting, 
having  been  used  both  subjectively  and  objectively  in  many  ways.  In 
numerical  weather  prediction,  climatological  statistics  are  used  in 
obtaining  the  weighting  functions  in  optimal  interpolation  objective 
analysis,  and  as  the  first  guess  in  some  iterative  objective  analysis 
techniques.  In  the  tropics,  climatology  constitutes  an  important  part 
of  statistical  forecasting  and  is  even  used  in  dynamical  prediction 
models.  In  very  short-range  terminal  forecasting,  climatological 
records  are  the  basis  of  the  widely  used  station  conditional  probability 
tables. 


Model  Output  Statistics  Technique.  Another  method  of  obtaining 
objective  forecasts  of  locaT weather  elements  involves  the  use  of  statis- 
tical models  to  complement  the  output  of  numerical  prediction  models. 

This  is  the  Model  Output  Statistics  (MOS)  technique  in  which  local  obser- 
vations of  weather  parameters  are  matched  with  output  parameters  from 
numerical  models  for  a period  of  one  year  or  more  [88],  Forecast  equa- 
tions are  then  derived  by  statistical  techniques  that  can  account  for 
biases  and  inaccuracies  in  the  numerical  model  and  local  climatology. 

The  MOS  technique  has  been  applied  to  forecasts  of  such  weather  elements 
as  precipitation,  clouds,  thunderstorms,  temperature,  and  wind  [89]. 


The  contributions  of  satellite  meteorology  to  modern  weather  prediction 
are  already  substantial.  Visible  and  infrared  imagery  data  and  radio- 
metric  atmospheric  sounding  data  from  satellites  are  currently  being 
applied  in  numerical  prediction  models.  Other  remote  sensing  tools, 
i.e.,  Doppler  radar  and  acoustic  sounders,  now  offer  interesting  new 
possibilities  for  improved  capability  in  the  prediction  of  severe  storms. 

Forecasting  of  Precipitation 

Precipitation  is  one  of  the  most  difficult  weather  elements  to  forecast, 
and  numerical  forecasts  of  precipitation  have  been  the  least  useful  of 
all  outputs  from  numerical  models.  Weather  forecasters  still  have  to 
depend  on  personal  skill  and  experience  when  making  precipitation  fore- 
casts, particularly  forecasts  of  precipitation  amount.  There  is  a 
seasonal  variation  of  skill  in  manual  quantitative  precipitation  fore- 
casts dependent  upon  the  predominant  scale  of  precipitation  systems, 
i.e.,  synoptic  scale  in  the  winter  to  the  small  mesoscale  in  the  summer. 
It  is  more  difficult  to  forecast  the  occurrence  and  distribution  of 
precipitation  with  thunderstorms  than  with  winter  storms. 
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Availability  of  Physical  and  Dynamical  Meteorological  Information. 

A great  deal  of  available  physical  and  dynamical  meteorological  infor- 
mation  has  application  in  precipitation  forecasting.  Calculations  of 
raindrop  growth  made  with  a parameterized  updraft  pattern,  has  shown 
that  when  the  maximum  drop  diameter  becomes  3 to  4 mm,  rainfall  can  be 
expected  to  grow  to  50  mm/hr  [90].  Calculations  of  water  content,  as  a 
function  of  the  temperature  of  the  cloud  atmosphere  and  mean  velocity  of 
ascending  air  in  the  cloud,  have  predicted  rain  shower  intensity  compu- 
tations with  a correlation  of  approximately  0.80  to  observed  intensity 
values  [91].  The  concentration  of  cloud  droplets,  i.e. , the  number  of 
nuclei  active  when  maximum  supersaturation  is  attained  in  a cloud,  is 
the  most  sensitive  factor  in  determining  rainfall,  and  there  is  a criti- 
cal drop  size  distribution  for  initiation  of  rain  [92]. 

Mesoscale  Effects  on  Precipitation.  Precipitation  intensity  usually 
fluctuates  with  periods  of  l0-20  minutes  and  2-6  hours,  which  are  con- 
sidered to  be  related  to  convective  activity  and  mesoscale  disturbances, 
respectively  [93].  In  many  aspects,  mesoscale  disturbances  have  gravity 
wave  characteristics;  therefore,  it  is  important  to  clarify  the  inter- 
action between  convective  activity  and  mesoscale  disturbances.  Conse- 
quently, a Portable  Automated  Mesonet  (PAM)  was  used  to  identify  surface 
mesoscale  circulations  influencing  thunderstorm  development  during  the 
National  Hail  Research  Experiment  '76  [94].  PAM  consists  of  a trailer- 
mounted  base  station  and  a network  of  remote  sampling  stations  which 
sample  surface  mesoscale  data  synchronously.  Data  (pressure,  tempera- 
ture, humidity,  rain,  wind  direction,  and  speed)  are  averaged  locally 
and  transmitted  digitally  via  a telemetry  link  to  the  base  station  where 
real-time  data  from  the  entire  network  are  displayed. 

Furthermore,  after  about  20  years  of  big  mesometeorology  planning,  but 
little  systematic  implementation,  it  appears  that  some  dramatic  advances 
are  imminent.  Plans  are  now  being  formulated  for  a special  severe  storm 
satellite,  STORMSAT,  for  mesoscale  coverage  in  the  early  1980's,  and  the 
Severe  Storms  and  Mesoscale  Experiments,  SESAME,  will  play  a large  role 
in  providing  scientific  and  observational  support  for  understanding  and 
modeling  mesoscale  phenomena  and  assessing  their  predictability  [95]. 

Statistical  Applications  to  Precipitation  Forecasting.  Rapidly 
developing  hydrological  research  activities  require  more  Knowledge  of 
statistically  analyzed  time  series  data  of  precipitation  [96].  There 
has  been  increasing  interest  toward  fitting  a finite  auto-regression  to 
the  time  series  data  and  estimating  the  variance  of  a time  series  by 
calculating  the  spectrum  from  the  estimated  auto-regression  coefficients 
and  the  one-step  predictor  error  variance  [97].  Digital  filters  have 
become  a prime  tool  for  analyzing  sampled  time  series,  but  the  main 
problem  in  filter  design  is  to  obtain  short  operators  with  fairly 
accurate  frequency  characteristics  [98]. 

The  statistical  distribution  of  rainfall  amounts  over  relatively  long 
periods  (months,  years)  have  received  considerable  attention  in  the 
past,  but  for  moderate  to  small  time  intervals  (weeks,  days,  hours) 
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satisfactory  distributional  r.ocels  are  re  difficult  to  find  [99]. 

The  synoptic-statistical  approach  has  been  used  to  forecast  floods  and 
droughts  in  the  tropics  fcr  periods  ranging  from  one  month  to  years. 
However,  J.  Narnia s,  a US  pioneer  in  extended  range  forecasting,  has 
expressed  a hope  that  more  dynari’icai  methods  will  be  used  in  the  future 
[lOOL  Properties  of  the  Korean  rainy  season  have  been  investigated 
statistically  from  the  frequency  of  occurrence  of  a bad  weather  (cli- 
matic) index  and  the  seasonal  variation  of  the  polar  front  [1011. 

A very  recent  example  of  the  use  of  climatological  data  as  applicable  to 
hydrometeorology  involved  an  investigation  of  heavy  rainstorms  in  the 
northern  Illinois  area  [1021.  Based  on  historical  data  of  rainfall  and 
associated  weather  conditions,  frequency  distributions  of  point  rainfall 
were  presented  for  rain  periods  ranging  from  5 minutes  to  72  hours  and 
recurrence  intervals  of  6 months  to  50  years.  Also  presented  were  the 
relationship  between  point  and  area  mean  rainfall  frequencies,  time  and 
space  distribution  characteristics  of  heavy  rains,  orientation  and 
movement  of  storms,  the  frequency  distribution  of  storm  centers,  synoptic 
weather  conditions  associated  with  severe  rainstorms,  and  the  diurnal, 
monthly,  and  seasonal  distribution  of  flood-producing  storms. 

A stochastic  model  has  been  developed  to  determine  the  cumulative  distri- 
bution function  of  the  total  rainfall  per  event  based  upon  a particular 
random  structure  of  space-time  rainfall  [103].  Two  random  variables  of 
the  spatial  rainfall,  i.e. , the  cumulative  rainfall  within  a season  and 
the  maximum  cumulative  rainfall  per  rainfall  event  within  a season,  were 
considered.  One-minute  and  4-minute  precipitation  rate  records  have 
also  been  used  to  develop  stochastic  model  relations  between  clock  hours 
and  1-minute  (or  4-minute)  precipitation  rate  distributions  for  use  in 
estimating  seasonal  and  annual  distributions  of  instantaneous  precipita- 
tion rates  for  locations  where  distributions  of  clock-hour  rates  are 
known. 

Variational  analysis  of  hourly  observations  has  led  to  recognition  of 
the  importance  of  the  planetary  boundary  layer  in  severe  storm  develop- 
ment. Vertical  momentum  transport,  frictional  convergence,  and  gravita- 
tional motion  have  been  found  to  be  associated  with  the  formation  of 
squall  lines  and  thunderstorms.  It  has  been  found  that  in  a barotropic 
atmosphere  where  convection  is  predominantly  surface-forced,  a high 
correlation  exists  between  surface  convergence  and  the  predictability  of 
rainfal 1 . 

US  National  Meteorological  Center  Boundary  Layer  Model.  The  NMC 
has  experimented  with  a Planetary  Boundary  Layer  (PBL)  model  (similar  to 
the  US  Air  Force  Global  Weather  Central  (GWC)  operational  model)  which 
has  eight  levels  between  the  surface  and  1600  meters  and  a horizontal 
grid  mesh  of  190.5  km  over  the  contiguous  US  [104].  Forecasts  based  on 
0000  GMT  data,  utilizing  all  conventional  data  sources  plus  satellite 
input,  were  run  to  2400  GMT;  results  showed  that  the  model  has  potential 
for  providing  usable  forecasts  of  rain  versus  snow  delineation,  cloud 
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ceiling,  and  areas  of  severe  weather  over  the  Eastern  two-thirds  of  the 
US.  However,  the  nodel  did  not  perform  well  during  the  summer  months 
and  over  the  mountainous  region. 

Numerical  quantitative  precipitation  forecasts  have  been  computed  as  a 
function  of  the  vertical  velocity  and  humidity  distribution  in  the 
atmosphere.  The  orographic  influence  on  the  vertical  velocity  and  the 
relation  between  vertical  velocity  and  static  stability  of  the  atmosphere, 
as  an  important  factor  in  the  condensation  process,  have  been  introduced 
into  the  forecast  model  [1051. 

The  release  of  isolated  sur.irier  showers  has  been  related  to  influences  of 
the  urban  heat  island  [106].  The  forecasting  of  precipitation  types, 
utilizing  the  1000-500  mb  thickness,  has  been  enhanced  significantly  by 
also  considering  the  mean  lapse  rate  in  that  layer  [107].  The  total 
energy  of  instability  in  the  700-500  mb  layer  plays  a fundamental  role 
in  the  rain  intensity  process  within  cumuliform  clouds  [108]. 

Another  technique  vjas  developed  for  utilizing  existing  series  of  satellite 
cloud  photographs  to  estimate  total  daily  precipitation  (rain  and  snow) 
over  a specific  drainage  basin  (Flathead  River,  Montana)  [109].  Basic 
data  consisted  of  surface  observations  of  daily  precipitation  over  the 
climatological  network  and  a series  of  one-per-day  photographs  (both  as 
large-scale  mosaics  and  individual  frames)  from  polar  orbiting  satellites. 
A high  correlation  coefficient  was  obtained  between  estimated  and  observed 
precipitation  amounts. 

Precipitation  Forecasting  in  Other  Nations 

a.  Soviet  Union 

There  has  been  much  activity  in  the  Soviet  Union  relating  to  preciDita- 
tion  forecasting.  One  operational  model  uses  the  equation  of  mot’m  and 
the  first  law  of  thermodynamics  to  derive  prognostic  equations  in  a 
baroclinic  atmosphere  applicable  to  the  short-range  (48  hours)  prediction 
of  isobaric  contour  patterns  and  vertical  velocity  from  which  precipita- 
tion amount  or  intensity  can  be  predicted.  In  another  approach,  the 
entry  time  for  condensation  and  the  quantity  of  combined  water  vapor  for 
a given  time  interval  are  determined  by  calculating  the  derivatives  of 
pressure  and  the  initial  values  of  temperature  and  moisture  at  different 
atmospheric  levels  [110].  Vertical  velocities  at  1 and  3 km  heights  are 
then  calculated  based  on  ground  pressure  charts.  The  quantity  of  con- 
densed v/ater  vapor  in  an  air  column  6 km  high  over  the  time  interval 
under  consideration  is  expressed  in  millimeters  of  precipitation  and 
used  to  forecast  precipitation  amounts  for  a period  of  12  hours. 

An  integral  circulation  index  which  numerically  evaluates  the  atmo- 
spheric circulation  of  the  Northern  Hemisphere  has  been  used  to  forecast 
the  formation  of  precipitation  over  Central  Asia  with  a lead  time  up  to 
6 months  [111],  Also,  an  equation  which  considers  the  influence  of 
solar  activity,  types  of  atmospheric  circulation,  and  a 2-year  periodic- 
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ity  in  Northern  Hemispheric  circulation  has  allowed  calculations  of 
monthly  and  seasonal  precipitation  amounts  over  Western  Asia  with  a 
claimed  reliability  of  about  90%  [112]. 

In  1975,  the  Soviet  Union  issued  accuracy  results  of  forecasting  monthly 
precipitation  totals  in  the  past  10  years  for  regions  of  Kazakhstan. 
Forecasting  accuracy  of  the  measured  monthly  precipitation  totals 
throughout  the  territory,  and  by  seasons  of  the  year,  was  distributed 
nonuniformly  (from  54  to  80%)  [113]. 

b.  Other  Techniques 

Elsewhere,  a numerical  model  which  improves  the  original  Bushby-Thompson 
10-surface  atmospheric  model  by  more  realistic  considerations  of  such 
processes  as  friction  between  the  atmosphere  and  the  earth's  surface  and 
the  vertical  transport  of  heat  and  moisture  to  simulate  clouds  has  been 
used  successfully  to  predict  precipitation  distribution  24  hours  in 
advance  in  the  British  Isles  and  Western  Europe  [114]. 

A Polish  technique  entails  the  use  of  statistical  methods  to  analyze 
satellite  TV  cloud  pictures  [115].  Every  cloud  system  on  TV  pictures  is 
described  by  a five-digit  word;  thus,  a definite  digit  in  a definite 
place  in  the  group  describes  a definite  feature  of  the  cloud  system 
appearance.  Using  this  method,  it  is  possible  to  correlate  cloud  system 
appearances  and  precipitation  phenomena;  regions  of  rain  and  showers  and 
their  intensity  changes  up  to  12  hours,  and  of  thunderstorms  up  to  24 
hours,  from  picture-taking  time  can  be  determined.  Cloud  system  appear- 
ances may  also  be  used  as  additional  information  in  quantitative  precipi- 
tation forecasts. 

Thunderstorm  Forecasting.  Since  its  introduction  to  meteorology, 
radar  has  been  used  extensively  to  observe  the  formation  and  movement  of 
thunderstorms.  It  is  sometimes  possible  to  accurately  predict  the 
passage  of  a widespread  squall  line  of  thunderstorms  several  hours  in 
advance  and  the  arrival  time  of  an  individual  thunderstorm  up  to  perhaps 
20  minutes  in  advance,  but  usually  these  times  are  practical  upper 
limits.  In  the  prediction  of  thunderstorm  paths,  the  possibility  of  the 
dissipation  of  existing  storms  or,  more  importantly,  the  development  of 
new  storms  in  the  area  ahead  of  existing  storms  must  be  considered 
[116]. 


a.  US  National  Weather  Service  Approach 

The  Technical  Development  Laboratory  (TDL),  NWS,  approach  to  the  short- 
range  forecasting  (0-2  hours)  of  thunderstorms  and  severe  weather  places 
dependence  upon  digitized  radar  data  to  depict  the  intensity,  movement, 
and  evolution  of  severe  weather  systems,  with  pattern  recognition  tech- 
niques being  used  to  isolate,  track,  and  predict  the  movement  of  signi- 
ficant radar  echoes  [117]. 
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TDL  has  also  developed  and  implemented  an  objective  method  of  2-6  hour 
rorecasting  of  the  probability  of  thunderstorm  occurrence  within  a 
square  area,  40-45  nmi  on  a side  Cl  18].  The  prediction  scheme  is  based 
on  a hybrid  combination  of  the  classical  statistical  and  MOS  approaches, 
with  the  probabilities  being  produced  by  multiple  linear  regression 
equations.  The  independent  variables  or  predictors,  e.g.,  wind  compo- 
nents, mixing  ratio,  potential  temperature  advection,  and  stability 
indices,  are  derived  from  routinely  observed  surface  atmospheric  variables, 
manually  digitized  radar  data,  localized  climatic  frequencies  of  thunder- 
storms, and  large-scale  numerical  model  forecasts  of  basic  upper-air 
variables.  Over  a short  forecasting  period  (spring),  forecasts  of  the 
climatic  (statistical)  frequency  and  persistence  of  thunderstorms  showed 
considerable  skill.  Subjective  examination  of  individual  cases  showed 
that  the  envelopes  of  high  probability  values  matched  well  with  the 
general  patterns  of  thunderstorm  occurrence. 

The  TDL  approach  to  medium  range  forecasting  (6-24  hours)  of  thunder- 
storms and  severe  weather  involves  the  combining  of  dynamics  and  statis- 
tics with  predictors  generated  from  meteorological  parameters  obtained 
from  NMC's  6-layer  PE  model  and  TDL's  own  3-dimensional  trajectory 
model.  Using  computed  trajectories  from  wind  forecasts  by  the  NMC 
model,  detailed  24-hour  forecasts  of  temperature  and  dew  point  are 
derived  by  computing  6-hour  variations  of  potential  temperature  and 
mixing  ratio  for  air  parcels  assumed  to  follow  paths  defined  by  the 
trajectories. 

The  TDL  trajectory  model  was  recently  modified  to  include  surface  land 
and  ship  reports  and  initialized  grid-point  values  from  the  NMC  six- 
layer  PE  model  [119].  In  the  new  analysis  procedure,  upper-air  and  PE 
data  are  assigned  relative  weights  by  means  of  a unique  function  that 
measures  the  asymmetry  of  the  radiosonde  station  distribution  at  the 
initial  position  for  each  trajectory. 

b.  Other  Methods  of  Forecasting  Thunderstorms 

Another  objective  technique  has  been  developed  for  modifying  12-  to  24- 
hour  precipitation  guidance  forecasts  from  the  NMC  PE  model  by  using 
manually  digitized  radar  and  the  LFM  model  12-  to  24-hour  precipitation 
forecasts  [120].  Developmental  data  were  from  two  stations  in  the 
Eastern  US  for  the  summer  season.  Constructed  radar  variables  were 
entered  into  a stepwise  multiple  regression  program  with  the  PE  precipi- 
tation probabilities  and  the  LFM  precipitation  forecast.  The  resulting 
equation  yielded  subjective  improvement  over  the  PE  guidance  amounts  by 
10-15%. 

In  the  objective  forecasting  of  the  probability  of  thunderstorm  occurrence, 
recognized  severe  storm  indicators  have  included  the  following:  (1) 
hydrostatic  quantities  (e.g.,  Showalter  Stability  Index,  moisture  con- 
tent, and  a combined  representation  ot  the  low-level  temperature-dew 
point  fields);  (2)  the  effect  of  differential  temperature  advection 
between  the  low  and  mid  troposphere;  (3)  kinematic  quantities  representing 
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the  presence  of  an  activating  mechanism  of  potential  Instability  (e.g., 
low  level  vorticity  genesis);  and  (4)  a destabilizing  distribution  of 
temperature  advectlon  between  low  and  mid  troposphere  with  the  simultan- 
eous presence  of  horizontal  convergence  In  low  levels  and  of  horizontal 
divergence  In  the  upper  levels  [1211. 

A Severe  Weather  Threat  (SWEAT)  Index  for  forecasting  severe  thunder- 
storms and  tornadoes  was  developed  at  6WC.  The  index  Is  empirical  and 
based  on  the  850  mb  dew  point,  sum  of  the  850  mb  temperature  and  dew 
point  minus  twice  the  500  mb  temperature,  speed  of  the  850  mb  and  500  mb 
winds,  and  the  directional  shear  between  the  850  mb  and  500  mb  winds 
[1221. 


c.  Sferics  Technique 

There  is  also  some  evidence  that  a technique  for  noting  the  occurrence 
of  lightning  discharges  and  their  bearing  from  a station  (sferics  sig- 
nals) can  be  used  for  thunderstorms  and  tornado  identification.  It  has 
been  observed  that  the  frequency  of  sferics  reaches  a pronounced  maximum 
prior  to  the  formation  of  severe  thunderstorms  or  tornadoes. 

A transistorized  lightning  counter  has  been  developed  in  the  People's 
Republic  of  China  and  used  for  forecasting  hailstorms.  Simple  In  struc- 
ture, easy  to  handle,  and  inexpensive,  the  instrument  is  suitable  for 
use  in  both  plains  and  mountain  areas  [123].  As  described,  when  a mass 
of  ha  11 -producing  clouds  appears  within  40  km  from  the  Instrument,  Its 
antenna  will  automatically  transmit  lightning-produced  signals  to  the 
counter.  A given  number  of  recorded  lightning  signals  Indicates  that  a 
hailstorm  is  impending. 

US  National  Hurricane  Center.  The  HURRAN  statistical  model  used  at 
the  US  National  Hurricane  Center  is  based  on  all  recorded  past  storms 
that:  (a)  occurred  within  a fixed  number  of  days  from  an  existing  storm, 
(b)  occurred  within  a certain  distance  from  the  storm  (usually  within  2.5 
latitude),  and  (c)  moved  with  similar  speed  and  direction.  This  permits 
a stratification  of  hurricane  climatology  for  use  in  forecasting  and  long 
range  planning  [124]. 

Other  purely  statistical  hurricane  prediction  models  used  by  the  US 
National  Hurricane  Center  (NHC)  include  CLIPER,  NHC  67,  and  NHC  72.  A 
statistical -dynamical  model  is  represented  by  NHC  73,  while  purely 
numerical  models  include  the  barotropic  SANBAR  and  baroclinic  Movable 
Fine-Mesh  (MFM)  model. 

Many  studies  of  the  joint  frequency  of  the  initial  and  final  conditions 
of  weather  elements,  e.g.,  rainfall  and  cloud  cover,  attest  to  the 
importance  of  the  initial  event  as  a predictor  of  the  later  event  [125]. 
By  assuming  the  Markov  process  (in  which  the  past  history  does  not 
affect  conditional  probabilities  of  events  defined  in  terms  of  the 
future)  and  a bivariate  normal  distribution,  a two-parameter  model  was 
developed  by  US  Air  Force  Researchers  which  charts  the  parameters  as 
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direct  functions  of  the  probability  of  the  Initial  event;  the  model  can 
then  be  used  to  estimate  conditional  probabilities  of  both  frequent  and 
rare  events. 

Previous  models  for  estimating  the  conditional  probability  of  an  event 
have  used  a categorized  event  as  the  Initial  condition,  e.g.,  "no  rain" 
or  “overcast"  at  zero  time.  However,  Initial  conditions  frequently  are 
known  In  detail  and  this  model  applies  those  detailed  conditions  In 
determining  conditional  probabilities  Cl 261. 

Probability  of  Precipitation  Forecasts.  Probability  of  Precipita- 
tion CPOP)  values  derived  from  estimates  of  the  gamma  probability  func- 
tion (as  developed  by  Thorn)  have  been  used  to  depict  regional  precipi- 
tation characteristics  over  the  Eastern  US,  and  Indicate  the  Importance 
of  lake  effects,  sea  breeze  convergence,  coastline  geography,  and  proxim- 
ity to  storm  paths  In  POP  predictions  C127]. 

Radar  Information  has  been  used  with  synoptic  and  radiosonde  data  to 
produce  more  accurate  short-term  POP  forecasts  Cl 28].  The  radar  provides 
details  about  area  coverage,  echo  shape,  number,  size,  orientation,  and 
movement  while  nonradar  parameter  predictors  Include  the  mean  relative 
humidity  In  the  convective  cloud  layer,  slice  stability  Index,  and  sub- 
cloud water  vapor  capacity. 

Regression  techniques  for  obtaining  POP  over  the  Western  US  from  500-mb 
flow  typing  have  represented  a substantial  Improvement  over  climatology 
when  the  flow  Is  significantly  different  than  climatological  C129].  In 
computing  correlations  between  meteorological  factors  and  snowfall,  most 
statistical  multiple  regression  equations  appear  as  a combination  of  the 
predictor  related  with  temperature  producing  heavy  snow  on  the  large 
scale  and  the  predictor  related  with  surface  wind  controlling  the  local 
distribution  of  snowfall  [130]. 

Elevation,  latitude,  date,  850-mb  temperature,  1000-500-mb  thickness, 
and  boundary  layer  temperature  predictors  have  been  used  to  produce 
objective  (computer)  POP  predictions  of  snow,  sleet,  or  freezing  rain  up 
to  48  hours  [131]. 

Distinguishing  variants  of  the  disposition  of  pressure  systems  producing 
winter  snowstorms,  together  with  data  of  snowstorm  frequency  In  specific 
regions  of  the  Arctic  during  particular  synoptic  processes,  have  been 
used  to  make  snowstorm  probability  forecasts  with  a lead  time  of  3-4 
days  [132]. 

Estimates  of  Probable  Maximum  Precipitation.  Estimates  of  Probable 
Maximum  Precipitation  (PMI*),  or  the  greatest  depth  of  precipitation  for 
a given  duration  with  no  allowance  made  for  long-term  climatic  trends, 
are  based  mainly  on  the  traditional  approach,  I.e. , moisture  maximiza- 
tion and  transposition  of  observed  storms  and  consideration  of  elevation, 
moisture  Inflow  barriers,  and  distance  from  the  moisture  source  [133]. 
Procedures  for  estimating  PMP  vary  with  amount  and  quality  of  available 
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data,  basin  size  and  location,  basic  and  regional  topography,  storm 
types,  and  climate. 

Forecasting  of  Cloud  Amount  and  Height 

US  Air  Weather  Service  l^^dels.  Numerical  methods  have  become  a 
popular  means  of  predicting  fields  of  layered  cloudiness.  One  of  the 
better  known  mathematical  formulations  in  cloud  forecasting  is  the  GWC 
model  consisting  of  three  modules,  i.e. , macroscale  clouds  (MSC),  5- 
layer  (5  LYR),  and  high  resolution  cloud  prog  (HRCP)  [1341.  Forecasts 
are  based  upon  three-dimensional  parcel  displacements  computed  from 
dynamic  numerical  prediction  models.  The  most  detailed  HRCP  produces  a 
3-,  6-,  or  9-hour  cloudiness  forecast  for  15  layers  from  the  surface  to 
55  Kft  at  a 25  nmi  grid  interval  covering  the  Northern  Hemisphere.  The 
other  models  produce  less  detailed  cloudiness  forecasts  for  periods  up 
to  48  hours. 

GWC  also  has  an  operationally  useful  convective  parameterization  tech- 
nique applicable  in  numerical  cloud  forecasting  [1353.  This  technique 
uses  certain  dynamic  variables  from  the  GWC  data  base  and  relates  these 
variables  to  statistical  convective  activity  through  a nomogram.  The 
amount  of  cumulus  cloud  cover  is  assumed  to  be  a function  of  stability 
and  the  subcloud  layer  moisture  convergence. 

Other  Techniques  for  Cloud  Prediction.  The  MOS  prediction  tech- 
nique  treats  cloud  ceiling  height  both  as  a categorized  and  as  a con- 
tinuous predictand  [1363.  Where  cloud  ceiling  height  is  categorized, 
the  Regression  Estimation  of  Even  Probabilities  (REEP)  screening  tech- 
nique is  used  to  develop  probability  forecast  equations.  Where  ceiling 
height  is  treated  as  a continuous  variable,  specific  ceiling  height 
forecast  equations  are  developed  by  ordinary  screening  regression. 
Forecasts  from  REEP  equations  are  generally  better  than  those  based  on 
specific  height  equations,  climatology,  and  persistence. 

The  numerical  six  hourly  output  from  the  6-layer  NMC  PE  model  has  been 
used  to  predict  average  daytime  cloudiness  in  spring  [1373.  Variables 
chosen  as  predictors  were  the  PE  model  predictions  of  mean  relative 
humidity  between  the  surface  and  500-mb,  the  700-mb  vertical  velocity, 
and  the  relative  humidity  trend  through  a 12-hour  period. 

An  objective  method  uses  satellite  data  from  the  visible  and  IR  window 
bands  to  derive  current  and  6-hour  forecasts  or  cloud  amounts  and 
ceilings  [1383.  Visible  data  were  carefully  calibrated  and  normalized 
for  solar  elevation  and  bidirectional  reflectance;  and  cloud  motions 
were  used  for  6-hour  forecasts.  Data  were  reduced  to  values  of  average 
brightness  and  its  standard  deviation  over  1/2  x 1/2°  areas  before 
analysis.  Results  have  indicated  a root  mean  square  error  for  current 
cloud  amounts  (on  a scale  of  0 to  1)  of  about  0.25  and  for  (6-hour 
forecasts  about  0.35.  The  occurrence  of  a ceiling  (more  than  0.5  cover- 
age) or  no  celling  was  specified  correctly  for  80%  and  70%  of  the  time 
at  0 and  6 hours,  respectively. 
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Qualitative  estimates  of  the  occurrence,  persistence,  and  dispersal  of 
low-level  cloudiness  over  Asia  were  obtained  by  applying  the  fundamental 
factors  of  convergence,  moisture,  surface  temperature  variation,  and 
nature  of  stratification  in  the  lower  1.5  km  layer  of  the  atmosphere. 

One-  to  two-hour  forecasts  of  the  variation  in  lower  cloud  boundary 
height  have  been  made  from  data  of  the  variation  rate  of  the  lower  cloud 
boundary  during  the  last  observation  period  as  determined  by  the  Laplacian 
transformation  of  pressure,  the  temperature,  and  relative  humidity  at 
the  ground  and  lower  cloud  boundary  height  [139]. 

Frequency  distributions  of  cloud  amounts,  as  described  by  simple  linear 
equations,  have  also  been  used  to  estimate  mean  values  of  cloud  cover. 

Fog  Forecasting.  Enhanced  IR  satellite  pictures  have  been  suggested 
as  a new  tool  for  forecasting  the  extent  of  warm  air  advection  fog  and 
stratus  formation  [1401.  Frequently  areas  in  which  fog  and  stratus  are 
most  likely  to  form  will  appear  as  relatively  dark  areas  on  the  satellite 
pictures  taken  a few  hours  after  sunset.  These  dark  areas  are  found 
normally  upwind  from  a moisture  source  and  appear  to  outline  the  boundary 
of  relatively  moist  air  in  lower  levels  of  the  atmosphere.  This  boundary 
can  be  monitored  and  used  as  an  estimate  of  the  extent  of  fog  and  stratus 
formation  during  the  night.  A method,  quick  enough  for  operational  use, 
has  also  been  developed  for  using  satellite-observed  cloud  brightness  to 
forecast  fog  and  stratus  dissipation  [141]. 

A two-dimensional  dynamical  model  has  been  developed  to  forecast  fog  for 
periods  of  up  to  8 hours  on  a local  scale  [142].  Time-dependent  winds 
are  generated  by  calculating  a stream  function  from  a model -produced 
vorticity  field,  and  a diurnal  surface  variation  is  simulated  by  sinu- 
soidally varying  the  lower  boundary  conditions  on  temperature.  Surface 
terrain  effects  can  be  incorporated  into  the  model;  however,  fine-scale 
data  are  required  and  wind  values  are  limited  by  an  inability  to  include 
synoptic-scale  variations. 

Another  fog  forecasting  method  is  based  on  the  dependence  of  the  rela- 
tive frequency  in  occurrence  of  visibilities  lower  than  1000  meters  upon 
the  daily  time,  surface  wind  direction  and  speed,  and  the  type  of  synoptic 
situation  [143]. 

A simple  method  of  forecasting  fog  at  Prague,  Czechoslovakia,  is  based 
upon  the  values  of  average  temperature  lapse  in  the  layer  between  the 
ground  and  850-mb  level  and  the  values  of  wind  speed  at  the  900-mb  level 
[144]. 

A conductive-radiative  model  has  been  used  to  predict  the  formation  and 
growth  of  radiation  fog  by  solving  numerically  the  heat  and  mass  trans- 
port equations  in  conjunction  with  an  approximate  form  of  the  radiative 
transfer  equation  [145]. 
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A technique  developed  in  the  Soviet  Union  for  predicting  fog  formation 
intensity  and  duration  is  based  upon  mathematical  determination  of 
expected  nocturnal  variation  of  the  specific  humidity  and  nocturnal 
minimum  temperature  in  the  atmospheric  surface  boundary  layer.  Verifi- 
cation was  claimed  in  88%  of  about  160  cases  in  which  fogs  were  noted 
during  the  night  and  weather  conditions  were  favorable  to  fog  formation. 

Forecasting  of  Temperature,  Humidity,  and  Wind  Speed  Profiles 

Temperature  Profile  Forecasts.  TDL,  NWS,  has  under  development  a 
large-scale  three-dimensional  planetary  boundary  layer  model  to  predict 
temperature,  humidity,  and  winds  in  the  lowest  2 km  of  the  atmosphere 
[146].  The  model,  consisting  of  12  levels  within  two  layers  (i.e.,  a 
50-meter  thick  surface  layer  having  constant  fluxes  of  heat,  momentum, 
and  moisture  and  a transition  layer  extending  to  2 km)  has  shown  good 
agreement  between  the  predicted  and  measured  boundary  layer  temperatures 
and  winds. 

The  Monin-Obukov  log-linear  model  has  been  examined  for  adequacy  in 
describing  temperature  and  wind  profiles  in  thermally  stratified  shear 
flows  and  diversified  thermal  stability  [147].  The  dimensionless  wind 
shear  and  lapse  rate  for  all  ranges  of  thermal  stability  examined  were 
shown  to  be  linearly  dependent  upon  the  dimensionless  height  derived 
from  the  log-linear  model. 

The  Galerkin  method  has  also  been  used  to  numerically  solve  the  non- 
linear initial  boundary  value  problem  describing  the  vertical  temperature 
profile  for  a thermally  coupled  soil-atmosphere  boundary  layer  in  a 
simple  physical  setting  [148].  However,  this  numerical  simulation 
requires  special  mathematical  techniques  for  handling  interior  boundary 
or  coupling  conditions. 

Humidity  Profile  Forecasts.  Data  from  11  widely  dispersed  geograph- 
ical locations  within  the  Us  have  been  sorted  with  respect  to  temperature 
intervals  of  5“C  and  relative  humidity  intervals  of  10%  for  the  lowest 
1400  meters  of  the  atmosphere.  These  data  have  been  assembled  in  a 
statistical  format  and  classified  into  altitude  Increments  of  200  meters 
for  each  site  by  season  of  the  year  [149]. 

Many  mathematical  descriptions  have  been  presented  for  profiles  of 
boundary  layer  air  temperature  and  wind  speed  over  various  types  of 
ground  vegetation  surrounded  by  irregular  and  inhomogeneous  terrain  [150]; 
while  the  humidity  profile  has  been  estimated  by  using  an  empirical 
formula  linking  changes  in  humidity  with  changes  in  temperature  and 
altitude  [151]. 

More  accurate  prognostic  stratification  curves  of  nighttime  temperature 
and  humidity  have  been  constructed  by  taking  into  account  the  nonlinear 
temperature  and  humidity  variations  in  the  surface  to  850-mb  layer  [1523. 
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Wind  Profile  Forecasts.  The  generalized  Ekman  equation  has  often 
been  used  for  micrometeorological  determinations  of  the  wind  profile 
[153].  Both  the  eddy  diffusivity  and  the  thermal  wind  are  important 
considerations  when  this  equation  is  applicable. 

It  has  been  claimed  that  wind  velocity  variations  with  height  can  be 
described  satisfactorily  by  the  expression  = Vq  (Z/Zq)'”,  where  and 

Vq  are  wind  velocities  at  heights  Z and  Z^,  respectively,  and  m * 0.125 
[1541. 

Comparisons  have  been  made  between  tower-observed  wind  profiles  and  pro- 
files predicted  from  models  based  on  mixing-length  theory  and,  in  the 
case  of  a single  roughness  change,  on  the  turbulent  energy  equation  [155]. 
Comparisons  with  the  mixing-length  model  are  moderately  good,  but  some 
features  of  the  observed  profiles  are  missing  in  the  theoretical  predic- 
tions. The  turbulent  energy  equation  model  has  given  slightly  less 
accurate  wind  profile  predictions. 

SUMMARY 

The  measurement  and  prediction  of  weather  variables  are  of  extreme 
importance  to  Army  decision  makers  in  tactical  planning  and  operations. 
Weather  changes  can  seriously  affect  weapons  systems,  target  acquisition, 
and  battlefield  surveillance  systems.  Weather  can  also  influence  the 
deployment  of  smoke  and  chemicals  and  selection  of  air  and  ground 
mobility  tactics. 


Measurement  Capabilities 

Remote  sensing  is  a relatively  new  tool  for  investigating  the  structure 
and  dynamics  of  the  atmosphere.  It  has  significant  applications  to 
weather  and  severe  storm  forecasting,  because  it  combines  reasonably 
accurate  data  with  spatial  and  temporal  resolution  (in  essentially  real 
time)  over  much  larger  volumes  of  the  atmosphere  than  is  possible  by 
other  means. 

With  the  arrival  of  the  meteorological  satellite,  a capability  for 
complete  global  weather  observations  was  realized.  Visible  and  infrared 
images  from  polar-orbiting  satellites  now  provide  day  and  night  sur- 
veillance of  weather  systems  over  the  entire  earth  [156].  Geostationary 
satellites  hovering  over  the  equator  now  send  data  back  to  earth  every 
half-hour  for  a substantial  portion  of  the  Western  Hemisphere  [157]. 
Moreover,  the  geostationary  satellite  has  proved  valuable  as  a means  of 
indirectly  providing  data  on  winds  in  remote  ocean  areas  through  its 
ability  to  track  cloud  motions. 

Radiometric  methods  have  been  developed  for  sounding  the  atmosphere  from 
satellites.  In  addition  to  the  IR  sounders,  which  measure  temperature 
and  moisture  profiles,  experimental  microwave  sounders  have  been  flown 
that  are  able  to  detect  precipitation  areas  within  large  cloud  masses. 
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Significant  improvements  have  occurred  in  some  areas  of  meteorology 
during  the  past  5 years.  One  particular  example  is  the  application 
of  vertical  temperature  profile  radiometer  data  to  large-scale  fore- 
casting at  the  National  Meteorological  Center  (NMC).  However,  routine 
applications  of  satellite  data  to  mesoscale  problems  are  still  confined 
largely  to  visual  and  IR  images  of  clouds,  although  quantitative  radiom- 
eter data  for  the  mesoscale  remain  on  the  verge  of  being  used  routinely. 

The  most  useful  ground-based  remote  techniques  involve  acoustics  and 
radar.  New  operational  uses  of  radars  have  not  increased  greatly  since 
1972,  although  Doppler  radar  techniques  are  progressing  rapidly.  Of  all 
the  surface-based  remote  techniques,  the  acoustic  radar  field  has  ad- 
vanced most  in  terms  of  routine  application. 

A very  good  ground-based  remote  sensing  capability  exists  for  character- 
izing wind  velocity  because  of  the  availability  of  a wide  range  of 
Doppler  techniques.  An  all-weather,  wind-profiling  capability  is  attain- 
able in  the  boundary  layer,  with  prospects  good  for  all-weather  profiling 
to  tropospheric  heights,  using  some  combination  of  microwave  and  IR 
Doppler  systems. 

The  ability  to  remotely  sense  temperature  and  humidity  is  much  poorer. 

It  seems  likely  that  it  will  be  practical  to  measure  temperature  and 
humidity  profiles,  perhaps  to  tropospheric  heights,  within  the  next  6 
years.  However,  it  does  not  seem  likely  that  remote  measurement  of  the 
three-dimensional  temperature  and  humidity  fields  will  be  available  from 
ground-based  remote  sensors  on  this  time  scale. 

The  development  of  weather  radar  following  World  War  II  brought  with  it 
a belief  that  it  would  produce  data  with  sufficient  time  and  space 
resolution  to  accurately  forecast  severe  thunderstorms.  But,  in  spite 
of  some  fine  technical  developments,  radar  has  not  made  the  impact  on 
severe  weather  forecasting  that  was  expected.  The  basic  problem  lies  in 
determining  how  severe  the  thunderstorm  is  at  any  given  time,  and  how 
severe  it  will  become  in  the  following  minutes  or  hours.  Present-day 
operational  radars  of  the  US  Weather  Services,  e.g.,  CPS-9,  WSR-57,  FPS-77, 
do  not  actually  measure  wind  speed,  hail  size,  or  rainfall  rate;  rather, 
wind,  rain,  and  hail  have  been  related  to  radar  return  intensity,  echo 
height,  and  echo  development  through  empirical  and  theoretical  studies. 

Recent  developments  in  electronics  and  data  processing  have  simplified 
the  task  of  obtaining  quantitative  data  for  estimating  these  meteorological 
parameters  through  the  use  of  digital  radar.  Digital  video  techniques 
offer  advantages  for  the  digital  processing  of  rainfall  data  for  flood 
and  trafficability  forecasting.  It  is  also  suited  for  producing  digi- 
tally Integrated  echo  Intensity  contours  of  severe  storms,  and  can  be 
conveniently  transmitted  over  land  lines  for  remote  display. 

A difficult  problem  still  remains  in  developing  relationships  between 
radar  measurements  and  intense  wind,  rain,  and  hail  due  to  the  scarcity 
of  meteorological  measurements  of  such  unusual  events.  Present  weather 
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radars  measure  intensity  and  tops  of  clouds  and  thereby  infer  storm 
severity,  whereas  radar  with  a Doppler  function  can  directly  measure  the 
phenomena  (wind)  that  comprise  the  threat.  Thus,  parent  circulations 
of  tornadoes  and  gust  fronts  can  be  directly  detected,  measured,  and 
predicted. 


Forecasting  Capabilities 

Two  technologies.  Numerical  Weather  Prediction  (NWP)  and  satellite 
meteorology  have  revolutionized  meteorology  in  the  past  20  years. 
Intermediate-range  forecasts  of  1-3  days  are  now  made  by  NMC,  while  the 
very  short-range  forecast  of  1-6  hours  is  becoming  more  impacted  by 
satellite  data  as  GOES  pictures  at  30-minute  intervals  become  available 
in  near  real  time.  Z'  merging  of  these  technologies  for  making  6-18  hour 
forecasts  through  the  use  of  mesoscale  hydrodynamic  numerical  models  is 
now  envisioned  with  other  tools  (e.g.,  mesocl imatology  studies  based 
upon  satellite  and  radar  data,  special-purpose  simple  models  or  forecast 
schemes,  and  empirical  relations  gained  from  forecaster  experience  and 
ingenuity)  seen  as  important  supplemental  but  important  forecast  aids. 

Current  numerical  prediction  model  capabilities  prove  that  significant 
progress  has  been  made  in  the  past  10-15  years  in  the  forecasting  of 
synoptic-scale  weather  features  and  weather  parameters  that  can  be 
derived  from  them,  but  the  progress  in  forecasting  precipitation  and 
other  small-scale  weather  elements  has  been  much  slower.  Fine-mesh 
numerical  models  have  resulted  in  in.proved  forecasts  of  the  location  of 
precipitation  areas,  but  important  problems  remain  to  be  solved  in  the 
application  of  numerical  techniques  to  the  forecasting  of  smaller-scale 
weather  phenomena,  such  as  thunderstorms  and  heavy  precipitation. 

Although  the  impact  of  operational  prediction  models,  from  the  barotropic 
to  the  primitive  equations,  has  been  most  noticeable  in  construction  of 
sea-level  and  500-mb  prognostic  charts,  progress  in  precipitation  and 
cloud  forecasting  based  on  numierical  weather  prediction  at  the  NMC  has 
been  modest  11581. 

Improvement  in  short-range  (6-24  hours)  forecasting  accuracy  is  expected 
to  continue  at  a deliberate,  if  not  slow,  rate  [159].  However,  it  is  not 
likely  that  POP  forecasts  will  undergo  any  dramatic  improvement  in  the 
foreseeable  future  [1601.  Also,  there  are  indications  that  the  effect 
of  truncation  and  other  errors  at  the  grid  points  of  numerical  models 
upon  the  exact  timing  of  large-scale  dynamics  is  very  serious.  This 
constitutes  a formidable  mathematical -physical  computational  problem 
which  must  be  solved  before  accurate  long-range  weather  predictions  can 
be  achieved  [161].  It  is  very  improbable  that  an  operational  forecasting 
ability  will  be  attained  within  the  foreseeable  future  for  the  prediction 
of  seasonal  mean  atmospheric  conditions  for  a year  in  advance  [1621. 

On  the  large-scale,  general  circulation  models  have  a role  in  studies  of 
climatic  processes  but  it  has  been  difficult  to  develop  models  whose 
climatologies  compare  favorably  with  observed  climatologies  in  any  great 
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detail.  Recognizing  this  problem,  the  World  Weather  Program  has  estab- 
lished a goal  that  would  provide  the  data  and  knowledge  needed  to  extend 
the  range  and  accuracy  of  weather  forecasts,  particularly  the  prediction 
of  climatic  fluctuations  and  changes.  This  action  is  very  timely  in 
view  of  expert  opinions,  based  upon  recent  climatological  studies,  which 
conclude  that  the  world  is  now  entering  a period  during  which  major 
climatic  change  is  apt  to  take  place  with  floods,  hurricanes,  droughts, 
tornadoes,  killing  frosts,  and  throttling  snowstorms  occurring  far  more 
often  and  randomly  during  the  coming  decades  than  ever  before. 
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